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Abstract. The extreme synchrotron BL Lac object H 2356-309, located at a redshift of z = 0.165, was observed from 
June to December 2004 with a total exposure of «;40h live-time with the H.E.S.S. (High Energy Stereoscopic System) 
array of atmospheric-Cherenkov telescopes (ACTs). Analysis of this data set yields, for the first time, a strong excess 
of 453 y-rays (10 standard deviations above background) from H 2356-309, corresponding to an observed integral flux 
above 200 GeV of I(>200GeV) = (4.1 ± 0.5) x 10~ 12 cirr 2 s~' (statistical error only). The differential energy spectrum of 
the source between 200 GeV and 1.3 TeV is well-described by a power law with a normalisation (at 1 TeV) of No = 
(3.00 ± 0.80 slal ± O.SLyJxlO^cnrV'TeV- 1 and a photon index of T = 3.09 ± 0.24 sla , ± 0.10 sys . H 2356-309 is one 
of the most distant BL Lac objects detected at very-high-energy y-rays so far. Results from simultaneous observations from 
ROTSE-III (optical), RXTE (X-rays) and NRT (radio) are also included and used together with the H.E.S.S. data to constrain 
a single-zone homogeneous synchrotron self-Compton (SSC) model. This model provides an adequate fit to the H.E.S.S. data 
when using a reasonable set of model parameters. 



Key words, gamma rays: observations - galaxies: active - galaxies: BL Lacertae objects: individual: H 2356-309 



1. Introduction 

The Spectral Energy Distribution (SED) of Active Galactic 
Nuclei (AGN) spans the complete electromagnetic spectrum 
from radio waves to very-high-energy (VHE; E > lOOGe V) y- 
rays. I n the widely-acce pted unified model of AGN (e.g. Rees 
1984; lUrrv & Padovanilll995l) . the "central engine" of these 
objects consists of a super-massive black hole (up to 1O 9 M ) 
surrounded by a thin accretion disk and a dust torus. In some 
radio-loud AGN, i.e. objects with a radio to B-band flux ratio 
Fsghz/Fb > 10, two relativistic plasma outflows (jets) presum- 
ably perpendicular to the plane of the accretion disk have been 
observed. 

AGN are known to be VHE y-ray emitters since the de- 
tecti on of Mrk42 1 above 300 GeV by the Whipple group 
JPunch et alj fl992). who pioneered the imaging atmospheric- 
Cherenkov technique. At very high energies, a number of AGN 
(^10) were subsequently detected by different groups using a 
similar technique. Almost all these objects are BL Lacertae (BL 
Lac) objects, belonging to the class of Blazars (BL Lac objects 
and Flat Spectrum Radio Quasars), i.e. AGN having their jet 
pointing at a small angle to the line of sight. The only confirmed 
VHE detection of an extragalactic object not belonging to the 
BL Lac class is the giant radio galaxy M 87 lAharonian et alJ 
l2003U2005dl) . 

Two broad peaks are present in the observed SED of AGN. 
The first peak is located in the radio, optical, and X-ray bands, 
the second peak is found at higher energies and can extend 
to the VHE band. The observed broad-band emission from 
AGN is commonly explained by two different model types. 
In leptonic models, the lower-energy peak is explained by 
synchrotron emission of relativistic electrons and the high- 
energy peak is assumed to result from inverse Compton (IC) 
scattering of electrons off a seed-photon population, see e.g. 
ISikora &' lVIadeiski d200lh and references therein. In hadronic 
models, the emission is assumed to be produce d via the inter- 
action s of relativistic pro tons with ma tter ( Pohl & Schlickeiser 
20Q0 |, ambient p hotons JMannheimll 19931) or magnetic fields 
( Aharonian 2000), or via the inter actions of relativistic pro tons 
with photons and magnetic fields (Miicke & Protheroe 2001). 

The observed y-ray emission from BL Lac objects shows 
high variability ranging from short bursts of sub-hour dura- 
tion to long-time activity of the order of months. Detailed 
studies of variability of BL Lac type objects can contribute 
to the understandi ng of their intrinsic acceleration mecha - 
msms (e.g. lKrawczvnskie 7al] l200lt lAharonian etall l2002). 
Additionally, observations of distant objects in the VHE 
band provide an indirect measurement of the SED of the 
Extrag a lactic Background Lig ht (EBL), see e.g. lStecker et alJ 
lll992l) : Prim ack et alJ d 19991) and references therein. Due to 
the absorption of VHE y-rays via e + e~ pair production with 
the photons of the EBL, the shape of the observed VHE spec- 
tra is distorted as compared to the intrinsically emitted spectra. 
Using a given spectral shape of the EBL, the observed AGN 
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spectrum can be corrected for this absorption. The resulting in- 
trinsic (i.e., corrected) spectrum can then be compared to basic 
model assumptions on the spectral shape of the y-ray emis- 
sion, thereby constraining the applied shape of the EBL. In this 
context, it is especially important to detect AGN at higher red- 
shifts but also to study the spectra of objects over a wide range 
of redshifts, in order to disentangle the effect of the EBL from 
the intrinsic spectral shape of the objects. To date, the redshifts 
of VHE emitting BL Lac objects with measured spectra range 
from z = 0.033 to z = 0.129. 

The high frequency peaked B L Lac object (HBL) 
H 2356-309, identified in the optical bv lSchwartz et alJJl989h . 
is host ed by an elli ptical galaxy located at a redshift of z — 
0.165 dFalomolll99ll) . The object was first detected in X-rays 
by the satellite experiment UHURU dForman et al.ll 19781) and 
subsequently by the Large Area Sky Survey experiment on- 
board the HEAO-I satellite JWood etalJll984. The spectrum 
of H 2 356-309 as observed by BeppoSAX Jcostamante et"all 
2001) is not compatible with a single power law model, indi- 
cating that the peak of the synchrotron emission lies within the 
energy range of BeppoSAX. A broken power law fit yields a 
synchrotron peak around 1 .8 keV, with a detection of the source 
up to 50 keV. These observations qualified the object as an ex- 
treme synchrotron blazar. 

A sel ection of TeV candidate BL La c objects was pro- 
posed by ICostamante & Ghisellinil J2002t) . The objects were 
selected from several BL Lac samples and using informa- 
tion in the radio, optical and X-ray bands. VHE predic- 
tions for the selected objects were g iven by the author s 
based on a parametrisation proposed by [Fossati et alJ {T998), 
suitable for predictions of high state flux of an average 
source population. The authors also gave VHE flux predic- 
tions based on a simple one-zone homogeneous SSC model 
jGhisellini et alJEo02l) . appropriate for a quiescent state of the 
specific VHE source candidate. H 2356-309 is included in this 
list and the predicted integral flux values above 300 GeV for 
H 2356-309 are 8.4xl0~ 12 cm _2 s _1 for the parametrisation and 
l^xlO^cirrV 1 for the SSC model. It should be noted that 
no absorption due to the EBL was taken into account in these 
calculations. 

In this paper the discovery of VHE y-rays from 
H 2356-309 with the H.E.S.S. Cherenkov telescopes in 2004 
is reported. With a redshift of z = 0.165, H 2356-309 is one 
of the most distant AGN detected at VHE energies so far. 
H 2356-309 was observed by H.E.S.S. from June to December 
2004 (see sections |2 and [3}- Simultaneous observations were 
carried out with RXTE (Rossi X-ray Timing Explorer) in X- 
rays on 11th of November 2004 (see section WA\ . with the 
Nancay decimetric radio telescope (NRT) between June and 
October 2004 (see section IO and with ROTSE-III (see sec- 
tion 14. 3> in the optical, covering the whole 2004 H.E.S.S. ob- 
servation campaign. 

2. H.E.S.S. Observations 

The system of four H.E.S.S. ACTs, located on the Khomas 
Highlands in Namibia (23°16'18"S 16°30'00"E), is f ully op- 
erational since December 2003. For a review see, e.g.. lHintonl 
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H.E.S.S. data are taken in runs with a typical duration 
of 28 minutes. The data on H 2356-309 were taken with the 
telescopes pointing with an offset of 0.5° relative to the object 
position (wobble mode, offset in either right ascension or dec- 
lination). The sign of the offset is alternated for successive runs 
to reduce systematic effects. H 2356-309 was observed with 
the complete stereoscopic system from June to December 2004 
for a total raw observation time of more than 80 hours. In or- 
der to reduce systematic effects that arise due to varying obser- 
vation conditions, quality selection criteria are applied before 
data analysis on a run-by-run basis. The criteria are based on 
the mean trigger rate (corrected for zenith-angle dependency), 
trigger rate stability, weather conditions and hardware status. 
During the 2004 observations of H 2356-309 atmospheric con- 
ditions were not optimal (due to brushfires), resulting in a dead- 
time corrected high-quality data set of *40 h live-time at an av- 
erage zenith angle of 20°. 

Th ese data are calibrated as described in lAharonian et all 
Thereafte r, before shower recons truction, a standard 
image cleaning JLemoine-Goumard et alJ l2005) is applied to 
the shower images to remove night-sky background noise. 
Moreover, in order to avoid systematic effects from shower 
images truncated by the camera edge, only images having 
a distance between their centre of gravity and the centre of 
the camera of less than 2° are used in the reconstruction. 
Furthermore, a minimum image amplitude (i.e., the sum of the 
intensities of all pixels being part of the image) is required 
for use in the analysis to assure a good quality reconstruc- 
tion. Previous H. E.S.S. publicati ons are mostly based on the 
standard analysis ( lAharonian et alJ2005al) . Here we present re- 
sults f rom the 3D Model analysis wh ich is presented in de- 
tail in | L emoine- Goumard et alJ J2005I) and was also used in 
lAharonian et alJ J2005bl) . This method uses independent cali- 
bration and simulation chains and is briefly described in the 
following paragraphs. 

2.1. H.E.S.S. 3D Model Analysis 

The principle of the 3D Model r econstru ction 

method JLemoine-Goumard & Degrangd [2004; 

lLemoine-GoumarTetTlTT2^0^riAliarorriarret alJ l2005bl) is 
based on a 3-dimensional (3D) shower model using the 
stereoscopic information from the telescopes. The shower 
is modelled as a 3D Gaussian photosphere with anisotropic 
angular distribution. For each camera pixel the expected light 
is calculated with a path integral along the line of sight. The 
observed images are then compared to the model images 
using a l og-likelihood fit with eig h t para meters, described in 
detail in ILemoine-Goumard et alJ d2005l) . For each detected 
shower at least two images are required for the reconstruction 
of the angle 6 (the angle between the object position and 
the reconstructed shower direction), shower core impact 
position (measured as a radius with respect to the center 
of the telescope array), energy and the transverse standard 
deviation cr T of the shower. The dimensionless reduced 
3D width a) (used for y-hadron separation) is defined as 
(jj = <Xtp/D s , with the density of air p and the column density 



Ds at shower maximum. The energy spectrum of the y-ray 
excess is obtained from a comparison of the reconstructed 
energy distributions to the expected distributions for a given 
spectral shape. For the determination of the expected number, 
y-ray acceptances calculated from simulations are taken into 
account. This forward folding me thod was first developed 
within the CAT collaboration (Piron. F. 2000). The quality of 
the reconstruction improves when only events with a number 
of triggered telescopes N te i > 3 are accepted. In the analysis 
used in this paper, results from the 3D Model are given using 
this cut. A spectral analysis using an N te i > 2 cut yields 
consistent results. 

The on-source data are taken from a circular region of ra- 
dius 9 around the object position (on-source region). The back- 
ground is estimated from 11 control regions (off-source re- 
gions) of the same size and located at the same radial distance 
to the camera centre as the on-source region and normalised ac- 
cordingly. The significance in standard deviations above back- 
ground (subsequently cr) of any excess is calculated following 
the likelihood method of iLi & Mai Il983h . All cuts applied are 
summarised in Tabled Cut values were optimised using simu- 
lated y-ray samples and independent background data samples. 



Table 1. Summary of applied cuts. The cut in image amplitude 
is given in photoelectrons (ph.e.). 



cut 


value 


distance 


2deg 


image amplitude 


>60ph.e. 


# of telescopes 


>3 


core impact position 


<300m 


mean reduced scaled width 




mean reduced scaled length 




reduced 3D width a> 


<0.002 


e 2 


<0.01 deg 2 



3. H.E.S.S. Results 

3.1. Signal 

In Figure^ the distribution of the squared angular distances 9 2 
from H 2356-309, reconstructed with the 3D-Model (N te i > 3), 
is shown. In the data taken from the on-source region a clear ac- 
cumulation of events is seen at low # 2 -values, i.e. close to the 
position of H 2356-309. The off-region shows a flat distribu- 
tion as expected for a pure background measurement. With a 
number of on-source events of N „ = 1706 off-source events 
N /f = 13784 and a normalisation factor o , =0.0909 (the ra- 
tio between the solid angles for on- and off-source measure- 
ments), the data yield an excess of N on - aN /f = 453 y-rays at a 
significance-level of 11.6 a. A fit of a 2-dimensional Gaussian 
to an uncorrected excess sky-map yields a point-like emission 
and a location (23/i59m09.42s± 2.89s, -30°37'22.7"± 34.5") 
consistent with the positio n of H 2356-3 09 (23/i59m07.8s, - 
30°37'38"), as obtained by Falomo ( 1991) using observations 
in the optical and near-infrared. These results are summarised 
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Fig. 1. Distribution of squared angular distances 9 1 for 
H 2356-309, using the 3D-Model reconstruction and N, e i > 3. 

in Table|2l Additionally, the results from the standard analysis 
are given. The results for the 3D Model analysis used in this 
paper are also given using N te i 2> 2 for easier comparison with 
the standard analysis. 



Table 2. Summary of analysis results from the 3D Model anal- 
ysis (using two different cuts in telescope multiplicity). For 
comparison, the standard analysis results are also given. 





3D Model 


standard analysis 




N tel > 3 


Nua > 2 


N td > 2 


N „ (events) 


1706 


4389 


3776 


N of f (events) 


13784 


40420 


35280 


normalisation a 


0.0909 


0.0909 


0.0903 


Excess (events) 


453 


715 


591 


Significance 


11.6(7 


10.9 a- 


9.7 cr 



3.2. Energy Spectrum and Variability 

The differential energy spectrum obtained from the 3D-Model 
analysis {Ntel ^ 3) is shown in Figure |2] The spectral pa- 
rameters were obtained from a maximum likelihood fit of a 
power law hypothesis dN/dE = No(E /TeV)~ r to the data, 
resulting in a flux-normalisation of No = (3.00 ± 0.80 sta t ± 
0.31 sys ) x 10" 13 cm _2 s _1 TeV _1 , and a spectral index of T = 
3.09 + 0.24 stat + 0.10 sys . The x 2 vame of the spectral fit is 
6.6 for 7 degrees of freedom, corresponding to a x 2 -probability 
of P(x 2 ) = 0.47. In order to eliminate any systematic effects 
that might arise from poor energy estimation at lower energies 
(due to an over-estimation, on average, of low energies), the 
beginning of the fit range is set to the value of the post-cuts 
spectral energy threshold, i.e. 200 GeV for the 3D Model anal- 
ysis of this data set. Systematic errors (0. 1 for the index and 
20 % for the flux) are dominated by atmospheric effects, i.e. a 
limited knowledge of the atmospheric profile needed as input 
for the simulations. A detail ed description o f systematic errors 
can be found in e.g., Aharonian et al. ( 2006). The parameters of 
the spectral fit are summarised in Table[5] Additionally, the re- 
sults from the standard analysis are given for comparison. The 
data-points used in Figure|2]are listed in Table[4] 




_, , , , , , , i \LJ 

0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Energy [TeV] 



Fig. 2. Differential energy spectrum of H 2356-309, as ob- 
tained from the 3D Model analysis. The shaded area gives the 
confidence region (1 standard deviation) for the spectral shape 
under the assumption of a power law. The upper limit (arrow) 
is given with a 99 % confidence level. 

Table 3. Summary of the parameters of power law fits 
(dN/dE = M)(£/TeVr r ). For comparison with the 3D Model 
analysis (1) used here, the results from the standard analysis (2) 
are also given. Both analyses yield consistent results. In addi- 
tion to the parameter values, the x 2 probabilities of the fits are 
given. 





r 






N 

[10- 13 cm-V^eV- 1 ] 


pix 1 ) 


(1) 

(2) 


3.09 
3.06 


± 0.24 stat 
± 0.21a* 


± 0.10 sys 
± 0.10 sys 


3.00 + 0.80 stat ± 0.31 sys 
3.07 ± 0.75 stat ± 0.31 sys 


0.47 
0.60 



Table 4. Differential flux for different energy bins. The upper 
limit is given for a confidence level of 99 %. 



E 


O 




AO 




[TeV] 


[cm * 


_ 

s 


'TeV-'] 




0.223 


2.06 x 10" 


n 


7.13x10" 


12 


0.270 


2.18 x 10" 


1 1 


4.36 x 10" 


12 


0.329 


1.04 x 10- 


1 1 


2.46 x 10" 


12 


0.403 


5.77 x 10- 


12 


1.48 x 10" 


12 


0.494 


2.12x10" 


12 


8.51 x 10" 


13 


0.604 


1.26x10" 


12 


5.45 x 10" 


13 


0.742 


5.50 x 10- 


13 


3.29 x 10" 


13 


0.912 


5.77 x 10- 


13 


2.31 x 10" 


13 


1.113 


<4.06x 10 


-13 







The average integral flux above 200 GeV in the year 2004 
(fitting with a fixed spectral index of 3.09) is /(> 200 GeV) = 
(4.1 + 0.5) x 10" 12 cm" 2 s-' (statistical error only). Light-curves 
of I(>200GeV) versus the modified Julian date (MJD) of the 
observation are shown in Figure[5]for two different time-scales. 
The monthly flux variation is shown in the upper panel and the 
average monthly flux from June to December of 2004 is shown 
in the lower panel. A fit of a constant yields no evidence for 
nightly variability (P(^f 2 )=0.18). 

The ASM (All Sky Monitor) shows no significant X-ray 
excess nor variability in the same monthly intervals. 
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-12, 



-1 




= 2.43 ± 0.11 and a flux of 9.7+°| x 10"" erg cm *s ' in 
The x 1 °f the fit is 11 for 14 degrees of 
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Fig. 3. Flux above 200 GeV from H 2356-309 as measured 
by H.E.S.S. (3D Model Analysis) in 2004. Upper panel: 
each point corresponds to the flux measurement in one night. 
Lower panel: averaged flux in six different time windows (June 
to December 2004). Both light-curves are consistent with a 
constant emission over time. The RXTE observations (see 
Section l4~Tl were carried out on MJD 53320. During that night, 
no H.E.S.S. observations could take place due to poor weather 
conditions. 



the 2-10 keV band 

freedom, or a^ 2 probability P(x 2 ) = 0.68. This flux level is 
approximately a factor of 3 l ower than what was repo rted from 
BeppoSAX observations bv lCostamante et aD J200ll) and with 
a softer photon index than what was observed above 2 keV. 



4.2. NRT analysis 

The Nan^ay radio-telescope is a single-dish antenna with a col- 
lecting area of 200 x 3 4.56 m 2 e quivalent to that of a 94 m- 
diameter parabolic dish Jvan Driel et alJl99^) . The half -power 
beam width at 11cm is 1.9arcmin (EW) xll.5arcmin (NS) 
(at zero declination), and the system temperature is about 45 K 
in both horizontal and vertical polarisations. The point source 
efficiency is 0.8 KJy , and the chosen filter bandwidth was 
12.5 MHz for each polarisation, split into two sub-bands of 
6.25 MHz each. Data were processed using the Nangay local 
software NAPS and SIR. 

A monitoring program with this telescope on extragalactic 
sources visible by both the NRT and H.E.S.S. is in place since 
2001. For the campaign described here it consisted of a mea- 
surement at 1 1 cm every two or three days. Between 4 and 14 
individual 1 -minute drift scans were performed for each ob- 
servation, and the flux calibration was done using a calibrated 
noise diode emission for each drift scan. The average flux for 
the measurements carried out between 1 1 June and 10 October 
2004 was 40 + 8 mJy. This observed flux is most likely domi- 
nated by emission produced in jet regions further out from the 
core and thus represents an upper limit of any emission model 
for the total SED. 



4. Multi-Wavelength Analysis and Results 

4.1. RXTE analysis 

The RXTE/PCA djahoda & PCA Team! Il996l) observed 
H 2356-309 twice for a total of 5.4 ks on 1 1 November 2004 
after a Target of opportunity was triggered on this target. Due 
to poor weather conditions, H.E.S.S. -observations were not 
possible on the night of November 1 1, or on the 2 prior nights. 
However, the RXTE observations can be considered simul- 
taneous to the H.E.S.S. observation campaign of November. 
The STANDARD2 data were extracted using the ftools in 
the HEASOFT 6.0 analysis software package provided by 
NASA/GSFC and filtered using the RXTE Guest Observer 
Facility recommended criteria. Only the signals from the top 
layer (X1L and X1R) are used from the PCA. The average 
spectrum shown in the SED (Figure |4]i is derived by using 
PCU0, PCU2 and PCU3 data. The faint-background model 
is used and only the 3-20 PHA channel range is kept in 
XSPEC v . 1 1 . 3 . 2, or approximately 2-10 keV. The column 
density is fixed to the Galactic value of Nh = 1.33 x 10 20 cm" 2 
obtained from the PIMMS nH program 1 and is used in an 
absorbed power law fit. This yields an X-ray photon index 

1 See http://legacy.gsfc.nasa.gov/Tools/w3pimms.html 



4.3. ROTSE-III Analysis 

The ROTSE-III (Robotic Optical Transient Search Experiment) 
array is a world-wide network of four 0.45 m robotic, auto- 
mated telescopes built for fast (« 6 s) response to GRB trig- 
gers from satellites such as HETE-2 (High Energy Transient 
Explorer 2) and Swift. The ROTSE-III telescopes have a 
wide (1.85° x 1.85°) field of view imaged onto a Marconi 
2048x2048 pixel back-illuminated thinned CCD and are 
operated without filters. The ROTSE-III systems are de- 
scribed in detail in lAkerlof etaD J2003h . The ROTSE-IIIc tele- 
scope, located at the H.E.S.S. site, has been used to perform 
an automated monitoring programme of bla zars, inclu ding 
H 2356 -309. Data is analysed as described in Ahar onian et alJ 
d2005ch and references therein. During the observation periods 
covered by H.E.S.S. the apparent R-band magnitude m(R) from 
H 2356-309 as measured by ROTSE-III has its maximum at 
m(R) = 16.1 and its minimum at m(R) = 16.9. The host galaxy 
has b een resolved in the optical dFalomolll99lt IScarpa et alJ 
2000) and near-infrared (Cheung et al J I2003I) . These observa- 
tions show that H 2356-309 is a normal elliptical galaxy with 
an effective radius of about 1.8 arcsec in the R band. The con- 
tribution of the galaxy to the observed ROTSE-III flux is esti- 
mated to be m(R) = 17 using a standard de Vaucouleurs radial 
profile. 
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Fig. 4. Spectral energy distribution of H 2356-309. Above 200 GeV, the results obtained in this paper are used. Data from NRT, 
ROTSE-III and RXTE are simultaneous to the H.E.S.S. observations and shown as filled symbols. All other data are archival and 
shown as open symbols. The ROTSE-III flux is given after galaxy subtraction. Radio emission arises from regions further out in 
the jet. A single-zone homogeneous SSC model described in the text is shown as a solid line. The dashed line shows the SSC 
model without absorption. The inlaid box shows a zoom in the VHE regime. 



5. Discussion 

In Figure 0] a broad-band SED obtained from archival data, 
and simultaneous optical (ROTSE-III) and X-ray (RXTE) data 
together with the H.E.S.S. results presented in this paper is 
shown. Additionally, the result of a simple leptonic model 
is given as a solid line. The simplest leptonic scenario is a 
one-zone homogeneous, time independent, synch rotron self- 
Comp ton (SSC) model as initially proposed by Ijones et all 
ill 974 for compact non-thermal extragalactic sources. Here, 
we adopt a description with a spherical emitting region of ra- 
dius R and homogeneous magnetic field B, propagating with 
Doppler factor 5 with respect to the observer. The high energy 
electron distribution is described by a broken power law be- 
tween Lorentz fact ors r m j n and ym„, with a break at yb and a 
normalisation K (IKatarzvriski et al.ll2.00ll) . This SSC scenario 
is used, taking into account absorption by the EBL, to repro- 
duce the simultaneous data of H 2356-309. The density of the 
EBL is not well known in the yum wavelength regime. Given 
the high redshift of H 2356-309 (z = 0.165) and a compara- 
tively hard VHE spectrum, important constraints on the EBL 
density can be derived f rom the H.E.S.S . data. T his question 
is addressed in detail in lAharonian et alJ J2005el) where con- 
straints on the density of the EBL are derived from H.E.S.S. 
observations of 1ES 1101-232 and H 2356-309. Here, we use 



the P0.45 parametrisation from this paper which is very close 
to the lower limit from galaxy counts. 

As shown in Figure |4] using a model with a reasonable 
set of parameters provides a satisfactory fit to the simultane- 
ous x-ray and VHE data. The emitting region is characterised 
by 5 = 18, B = 0.16G and R = 3.4 x 10 I5 cm. The electron 
power-law distribution is described by K = 1.2 X 10 cm -3 , 
Tmin = 10 3 , 7m ax = 3 x 10 6 . The Lorentz factor at break en- 
ergy yb is located at y^ = 2.5 x 10 5 to place the peak emis- 
sion in between optical and X-rays while providing a good 
fit to the H.E.S.S. data. We take the canonical index a\-2 for 
the low-energy end and found CK2=4.0 for the high-energy end 
so as to fit the observed X-ray power law spectrum. Lowering 
Tmin extends the fit to lower frequencies and enhances IC emis- 
sion in the Me V-GeV domain. Synchrotron self-absorption cuts 
off emission below IR frequencies when using low values of 
ymin. Radio emission arises from regi ons further out in the jet . 
Similar to the case of PKS 2155-304 (lAharonian et al.ll200.5ch 
we cannot exclude a possible contribution of such an extended 
region to the optical flux measured by ROTSE-III. This may 
soften some of the above-mentioned constraints. 

Although our VHE observations provide strong constraints 
on the physical parameters of single-zone SSC models, there 
is still some freedom of choice for the parameters that could 
be constrained further by a better understanding of the origin 
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of the optical emission, a better spectral coverage in the X-ray 
and sub-TeV region and the observation of possible variability.. 

6. Conclusions 

The high frequency peaked BL Lac object H 2356-309, lo- 
cated at a redshift of z — 0.165, was discovered in the VHE 
regime by the H.E.S.S. Cherenkov telescopes. Two different 
reconstruction and analysis methods were applied to the data 
both yielding consistent results. No strong evidence for vari- 
ability in the VHE band is found within the H.E.S.S. observa- 
tions. The same holds true in the X-ray band, where the object 
does not show any strong flux variability, neither in the ASM 
nor in the pointed observations. Additionally, the RXTE flux, 
observed simultaneously to the H.E.S.S. observations, is lower 
than the previously-measured BeppoSAX flux. This might in- 
dicate that our observations took place during a relatively low 
state of emission. 

For the first time, an SED comprising simultaneous radio, 
optical, X-ray and VHE measurements was made. A simple 
one-z one SSC model, taking into account absorption by the 
EBL ( Aharon ian et alJEo05el) . provides a satisfactory descrip- 
tion of these data. 

Given the high redshift of the object, the observed H.E.S.S. 
spectrum provides strong constraints on the density of the EBL 
jAharonian et alj E005e). Future observations of H 2356-309 
with H.E.S.S. will improve the accuracy of the spectral mea- 
surement and might also allow an extension of the observed 
spectrum to higher energies. This will provide further con- 
straints on the absorption of y-rays by the EBL. 
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